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ABSTRACT: The performance of parametric conversion in achieving number amplification 
and duplication is analyzed. It is shown that the effective maximum gains G. remain well 
below their integer ideal values, even for large signals. Correspondingly one has output 
Fano factors F. which are increasing functions of the input photon number. In the inverse 
(deamplifier/recombiner) operating mode, on the contrary, quasi ideal gains G. and small 
factors F, ~ 10% are obtained. Output noise and nonideal gains are ascribed to spontaneous 
parametric emission. 


1. INTRODUCTION 

The ultimate transparency of optical networks is essentially quantum-limited and any improvement 
beyond the standard performance depends on availability of nonstandard high quality quantum 
amplifiers. The photon number amplifier (PNA) and the photon number duplicator (PND) are the 
quantum devices which are needed in direct detection. 1 The PNA ideally should affect the state 
transformation 

\n) — |Gn> (1) 

for integer gains G and input eigenstates |n) of the number. Similarly, the PND, instead of 
amplifying the photon number, produces two copies of the same input state for eigenstates of 
the number, namely 

I n) — ► | n,n) * (2) 

Both devices are particularly suited to local area network environments, where the minimum loss for 
user-derivation is 3dB (in average), and transparency rapidly degrades with the increasing number 
of users. In such situation the PNA represents the ideal preamplifier to be inserted before each 
derivation, whereas the PND — which ideally realizes the quantum nondemolition measurement of 
the number — could itself be used as an ideal lossless optical tap. 

The PNA and PND could also be profitably used in the inverse operating mode, namely the 
PNA as a number deamplifier and the PND as a number recombiner. The number deamplifier could 
be used as a number squeezer , allowing production of subpoissonian states from coherent light; the 
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number recombiner, on the other hand, could produce novel nonclassical radiation from input twin* 
beams [an example of such application in production of phase-coherent states 2 is proposed in Ref. 

[3]]- 


The concrete realization of high quality PNA and PND for practical applications is an arduous 
task. As explained in Ref. [4], number conversion, in a way similar to the customary conversion, 
requires a medium with a x^ 2 * or X^ susceptibility, but here with a phase* dependent polarizability. 
More precisely, almost ideal number conversion can be achieved upon modulating the nonlinear 
susceptibility at a (<? — l)-submultiple of the wavelenght of the amplified mode, G being the integer 
gain [feasibility studies of number conversion using multiple quantum wells heterostructures are 
currently in progress®]. The required phase- dependent polarizability in a x^ or x* 3 * * medium may 
also be regarded as an intensity-dependent coupling for a x* G * or X^ G+l * medium (simply from polar 
decomposition of the boson field operators). This suggests that a gain-two PNA should be simpler 
to realize than a generic G > 2 amplifier. However, as also explained in this paper, the intensity 
dependent coupling should follow the power low (a*a + l)* 1 ^ 2 , a*a being the number operator of the 
amplified mode. Such a decreasing factor is essentially the (1 + J)~ l/2 saturating behaviour of a two 
level system effective susceptibility in the inhomogeneous- broadening limit, 5 but it is not obvious 
that this power low — which is obtained in a semiclassical context — could survive in the quantization 
procedure. 


The previous observations quite naturally lead to ask if the conventional conversion could 
somehow simulate the number conversion, and what would be the range of physical parameters 
where ideal behaviour is better approximated: this is the subject of the present paper. Quite 
unexpectedly (see for example Ref. [1]) we find that ideal behaviour is never approached, even in 
the limit of large input signals. The most striking result is that conversion is never complete and, 
therefore, the effective maximum gains G. remain well below their integer ideal values, even for large 
input photon numbers: quantum mechanics thus reveals its subtle nature even for large quantum 
numbers, here in form of noise in amplifiers [for a discussion on applicability of the correspondence 
principle in a different context, see Ref. [7]]. 

The inverse devices — namely the number deamplifier and the number recombiner — are better 
approximated by parametric conversion than the direct ones. We will show that ideal gains are 
achieved in the large-n limit, whereas Fano factors F. remain nonvanishing but small ( F m ~ 10%). 
Therefore, it seems that at present the devices which are simplest to realize concretely should be 
the number deamplifier and the number recombiner (even though probably the limited output noise 
of the deamplifier could not be satisfactory for applications as number squeezer). 

After presenting the theory of the ideal devices in Sect.2, the connections between the 
conventional and the number conversions are explained in Sect.3, where a simple mean field approach 
for analytical evaluation of the conversion time is also given. In Sect. 4 the announced numerical 
results on conversion times, effective gains and Fano factors are presented. In Sect. 5 we conclude 
with some remarks on the physical interpretation of the nonideal behaviour in terms of spontaneous 
parametric emission. 


2. THE IDEAL NUMBER AMPLIFIER/DUPLICATOR 

In the Heisenberg picture the ideal PNA corresponds to multiplication of the number operator by 
the integer gain G 

a* a —* Ga*a , (3) 

a being the annihilator of the amplified mode of the field. Because of the integer nature of a f a, 


72 



the deamplifier does not trivially correspond to replace G into G -1 in Eq. (3). Actually the ideal 
deamplification is the following 

® f ® -♦ [G^VaJ , (4) 

where [x] denotes the integer part of x. Asa consequence, even in the ideal case, the deamplification 
has an input-dependent effective gain G. 


G. = 


IG-'n) 


< G~ l 


(5) 


for n input photons, and G. G -1 for large n. [As an example, the case G = 2 is depicted in Fig. 
4.] In terms of the shift operator e+ : e+|n) = |n + 1), the transformation (4) is obtained as follows 


U - (i+f , 

where now (e+) G |n) = |n + G). In fact, the map (6) corresponds to the following 3 

flt - ®(G) > 

where is a boson operator creating G photons at a time, 8 namely 

tt (G)l n > * + !|» + 0) , 


( 6 ) 

(7) 

( 8 ) 


[ a (G)*a|c)] — 1 , [o(G) , J = Go(c) 


The explicit form of a| G) is 


(9) 


.t _ 

®(C) - 

and from Eq. (10) it follows that 


{ 


[G^nKn-G)! 

n! 



(10) 


®(G) a (G) = [^ _l « ,a ] . (11) 

which is the deamplification (4). The direct amplification (3) corresponds to the inverse 
transformation 


®(G) “» ® f (12) 

[see Ref. [3] for more details about these maps]. The transformations (7) and (12) are essentially 
permutations of two different types of boson. For commuting modes [a, c] = [a, e*] = 0 the permuting 
map a «-► c is realized by the Heisenberg evolution 


where 


PaP = c, PcP — a , 


P = P' 



(13) 

(14) 
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However, u <z uid ®(G) do not coinxnute, it is convenient to consider & simultaneous chsn^e of the 
field mode (namely the amplifier also converts the frequency or changes the field polarization). In 
this case the amplifying operator is 


P(C) = exp 



exp 



®(G) C + C to (C) 



(15) 


The operator (15),now attains the transformations 

%)(1 ® e+)P( C ) = (e+) G ® 1 . P(G) [(e+) G ® i] P(G) = 1 ® e+ » ( 16 ) 


where in the tensor notation 0\ ® O 2 the first entry is for the d mode and the second entry for the 
c mode. The Schrodinger evolutions of the number eigenstates corresponding to the amplifying and 
deamplifying operating modes are 

^(C)|0,ti) = |Gn,0) , (17) 

P iG )\n,0) = \G(G-U)AG- l n)) y (18) 

where (x) = x — [x] denotes the fractional part of x and |n,m) = |n) a ®|m) c . If one would consider 
only one mode in the above transformations — say a — a frequency conversion P( 1 ) is needed. In this 
case the evolutions (17-18) rewrite 

-P( 0 ) ^(1)1^7 0) = K?n,0> (1®) 

[G ! n]) , (20) 

whereas totally ignoring the mode c corresponds to trace the transformations (19-20) over this mode, 
adopting a density matrix representation for states. In this way nonunitary transformations for the 
reduced density matrix of the signal mode a are obtained, which do not preserve the Newmann- 
Shannon entropy: these are the ‘photon fractioning’ and ‘multiphoton ’ transformations of Refs. [3, 9]. 
The mode c is responsible of the added noise which is present even in the ideal case (see Eq. (5)) 
and corresponds to the ‘idler mode’ of the customary linear amplification. 3 

Apart from the * phase shift — which can be obtained by changing the optical path of the 
b mode and which, however, for an input number eigenstate corresponds to an irrelevant overall 
phase factor — the evolution operator (15) comes from the interaction Hamiltonian in the Dirac 
picture 

Hf = a| c) c -I- h.c. (21) 


for a dimensionless evolution time 


T 

Tm ~ 2 


The Hamiltonian (21) has the following constants of motion 

j .4 = a f a + Gc'c , 

d A = G{G~ l a'a) = G{G~ l i A ) , 


( 22 ) 

(23) 

(24) 


and, because of identity (24), only s A must be specified. In the following evaluations we use the 
basis of the Hilbert subspace corresponding to fixed i. 4 eigenvalues 

| n), A = |s .4 - Gn,n) , n = 0, 1, ...[G _1 s,tl (25) 
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In this basis Eqs. (17-18) rewrite 


•P(G)|rc)Gn = |0)c„ , (26) 

%)|0)„ = |[C?- 1 n]) B , (27) 

whereas the Hamiltonian (21) takes the tridiagonal form 

Hi\n), = aj^ln- 1), + ajfl il» + 1)» , (28) 

7»( [C?~ 1 « J - n + 1) . (29) 

Conservation of the interaction Hamiltonian (21) itself corresponds to the resonance condition 
GijJ a — w c . In the nonresonating case a third pump mode d is needed with u = Gu? a - a/ c : Eq. (21) 
is obtained from the interaction Hamiltonian in the Schrodinger picture 

H f = a| G) cd + h.c. (30) 


in the parametric approximation of classical undepleted pump, namely with d in a highly excited 
coherent state. 

The photon number duplicator in some respect is similar to the gain-2 photon number amplifier. 
Instead of amplifying the number of photons, it produces two copies of the same input state for 
eigenstates of the number operator. If the input copies are carried by the modes a and b whereas 
the output by c, the duplication map reads 


|0,0,n) -► |n,n,0) (31) 

and is trivially inverted for n a = n*> [for the general case see Ref. [3]]. The state transformation 
(31) corresponds to the Heisenberg evolution 

e+ 9 e+ ® 1 i ® i ® , (32) 

which is obtained as permutation of the boson operators a( M ) and c, where now denotes the 
two-mode creator 

fl (U)l n «* w *) = \/(nun{n a ,7ifc}) + l|n a + l,n 6 + 1) , 

[«(i.o> afi.i)] = i , i«(i.i)i a<o + = 2a (i.i) • 

The foUowing realization of a|, ^ is obtained in Ref. [3] 

= ttV — rains-i — iiii , 

y/maxfo^Oji^} + 1 

In a way analogous to the PNA, the Dirac picture interaction Hamiltonian of the PND is 

Hi = a|j 1} c + h.c. , (36) 

with constants of motion 

io = ^ (a*a + b'b + 2c f c) , (37) 

d D = a f a - 6*6 . (38) 


(33) 

(34) 

(35) 
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The Hilbert subspace of interest for duplication corresponds to dp = 0; the subspaces for fixed 
eigenvalues sp are spanned by the eigenvectors 

|n), D = \»d -n,s D -n,n) . (39) 

For fixed s D the Hamiltonian (36) has the tridiagonal form 

H,\n.) =/#> |n -1), + Pn+\ l n + 1). , ( 4 °) 

0$ = \j{* - n + 1)» . (41) 

Frequency conversion and simultaneous duplication require a classical undepleted pump mode d at 
frequency = u> Q + with interaction Hamiltonian 

H f = + h.c. (42) 

3. NUMBER-OPTIMIZED DOWNCONVERSION 

The Hamiltonians (30) and (42) are complicated by the occurrence of the multiboson operators a' {C) 
and a|, An outlook at Eqs. (10) and (33) reveals that the <7-photon amplification corresponds to 
a susceptibility and the duplication to a x^* the followings the G = 2 case the simplest 

to attain in practice — will be considered only. For (a* a) ^ 2 the two photon operator a^ can be 
approximated as follows 

a| 2) ~ a t2 [2(a f a + 1)] * , ((<**«) > 2) • (43) 

On the other hand, for dp = 0 the two-mode operator oj, is simply 

tt | w) -«V« + ir*» t . (*'a = b'b). (44) 

Hence the Hamiltonians (30) and (42) become 

H'~a n [2(a t « + l)]"*ai + h. C . (PNA) , (45) 

H' = a'{a'a + l)-h'cd + \xx. (PND). (46) 

As a crude approximation we substitute the intensity-dependent factors in Eqs. (45-46) with their 
constant average values and use the customary four wave mixing Hamiltonians 

Hfivm = a t2 cd + h.c. (PNA), (47) 

HnvM = a'b'cd + h.c. (PND) . (48) 

In the parametric approximation of undepleted classical pump <f, Eqs. (47-48) correspond to the 
interaction Hamiltonians 

Hi = a ,2 e + h.c. (PNA), (49) 

H, = a f 6*c + h.c. (PND) . (50) 
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The interaction time is rescaled by y/Tj, Id being the intensity of the pump d: the relation between 
the dimensionless time r and the real time t (namely the length of the nonlinear medium) now reads 

T = X {3) \fidt • (51) 

The Hamiltonians (47) and (48) were already suggested by Yuen 1 , who inferred the amplifying 
performance from the conservation laws (23-24) and (37-38), with the assumption of complete 
conversion of the input signal. However, we will show that complete conversion is never achieved, 
apart from the case of one input photon. As an example in Fig. 1 the average output photon 
number is plotted for Hamiltonian (48) versus the interaction time r, for both cases: number 
duplicator ((n o ) 0 = (th,) 0 = 0, (n c ) 0 = n,) and number recombiner ({» a ) 0 = {ns) 0 = n,,{n c ) = 0). 
An oscillatory quasiperiodic (or long-time periodic) behavior is evident, conversion never being 
complete at any time: the ideal gain is not reached, and the unconverted photons contribute to 
the output noise. Therefore, the saturating factors in Eqs. (45-46) are crucial to get complete 
conversion. Semiclassically a similar saturating behavior oc (1 + J)~i is obtained for interaction 
of radiation with a two level system in the inhomogeneously-broadening limit or in the adiabatic- 
following regime: 5 however, a full quantum treatment is still lacking and would require a wideband 
analysis. Here we only consider the performance of parametric Hamiltonians (49-50) in achieving 
approximate PNA and PND. In this case the interaction time r = t. for conversion depends on the 
input photon number n, 


71 =/ ( ” c) ° ’ 

' 1 («.)o . 


({«o)o = 0 : direct operating mode) , 
((n c )o = 0 : inverse mode) , 


which, in order to simulate the intensity-saturating low in Eqs. (45-46), should behave as follows 


The conversion time (53) could be obtained tuning the pump intensity on the input photon number 
n,: for n, varying in a wide range, this would require a suited feedback mechanism based on a 
quantum nondemolition measurement of »,. In the following we give more accurate evaluations of 
t„, using either analytical methods (a mean field approximation) and numerical calculations. The 
results obtained in the two ways will be compared and discussed in the end. 


3.1 A mean field approximation 

In Ref. [10] a linearization procedure for parametric conversion has been proposed, where 
Hamiltonians (49-50) are approximated in a selfconsistent way by the ideal ones (21) and (36). 
As we will see in the followings, this approach is correct only in the limit of large input photons 
numbers in the amplified/duplicated channels (i.e. a and b modes), namely it is suited to describe 
the inverse operating mode only. The method allows evaluation of the conversion time t.: its major 
limitation is that it leads to exact conversion and, therefore, there is no systematic way to estimate 
quantum fluctuations and nonideal gains. As a consequence, the direct operating mode cannot be 
described in terms of the time-reversed transformation of the inverse mode, because in this case 
knowledge of the output noise become essential. Therefore, in this section we analyze only the 
deamplifier/ recombiner case. 

The starting point of the method is to rewrite Hamiltonians (49-50) in a form similar to the 
ideal ones (21) and (36), namely 

Hi = f(a'a)Ac + A i f(a i a)c , (54) 
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where 


and 


A = 


{ 


0 ( 2 ) 

“( 1 . 1 ) 


(PNA) , 
(PND) 


f (2. + 3 -(-I)*)* (PNA), 

nx) ~ { (x + 1)* (PND). 


(55) 


(56) 


The operator function /(a f a) will be treated as a c-number time-dependent effective coupling, to 
be determined selfconsistently a posteriori. The Hamiltonian (54) is rewritten as 

Hi = /(a^Ac* + /(a f a - v)A i c , (57) 

where, in order to have a unified description of the two devices, the integer number v is used 

(PNA) , (58) 

\2 (PND) . V ' 

We write a mean field Hamiltonian taking the intermediate value /( o*a) = /( a f a - * ) between the 
two forms in Eq. (57) and averaging on the input state. One obtains 

H„ F = f(n a (r))[Ac' + h.c] , (59) 

where 

n a (r) = (o t (r)a(r)) 0 , (60) 


and 



(61) 


(the oscillating (-1)"‘ (t) term in Eq. (56) is neglected). In the Dirac picture the time evolution of 


an operator O is written as follows 

0{t) ~ exp ( iH av r ) O exp {-iH av rj , (62) 

using the time-averaged Hamiltonian H av 

H av = ^ (Ac* + A'c) , (63) 

6(t)= f T /(Tl 0 (r , ))dr , . (64) 

Jo 

The evolution of the operators A and c takes the simple form 

A(r) = Acos0(r) + tcsintf(r) (65) 

c(t) = e cos 0(r) + iA sin 0(r) . (66) 
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Figure 1: Time evolution of the output signal (n„) (figures on the top) and of the r.m.s. output 
noise (Ah 2 ) (figures on the bottom) for parametric conversion (Hamiltonian (50)) of input number 
states with n, = 10. The two figures on the left refer to the number duplicator ((n o ) 0 = (nfc)u = 
0, (n c ) 0 = n,); those on the right to the number recombiner ((»„)o = ( «6>o = n,,(n c ) = 0). The 
small circles enclose the conversion point corresponding to r = t.. 


We are now in a position to evaluate n a (r) selfconsistently. From Eq. (65) one has 

n^i(r) = (i4 t (r)y4(r))o = (A'A) 0 cos 2 9(t) . 

For large input photons n, and v = 1 the expectation n.-i can be approximated as 

n.4 = K/2] - n a /2 . 

From Eqs. (64) and (67) we obtain the following integral equation for n„(r) 

n a (r) = n, cos 2 6(t) = n, cos 2 /(n a (r'))dr' . 

Differentiation of Eq. (69) leads to 

s/vdr = d6 ^n, cos 2 9 + ^ 

From Eqs. (65) and (69) one can see that complete conversion occurs at r = r. such that 

Hr.) = l . 


(67) 

( 68 ) 

(69) 

(70) 

(71) 
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Figure 2: The beat conversion time r m for parametric Hamiltonians (49-50) (PNA on the left 
and PND on the right). The squares are for the amplifier /duplicator, the circles for the 
deampliiier/ recombiner. The lines without dots represent the mean-field approximation. 

After integrating Eq. (70) from 9 = 0 to 9 = * we find the conversion time as a function of the 
input photon number n t 




where K(k) denotes the complete Jacobian elliptic integral 


w 

K(k) = / ( 1 - fcsin 2 ®)“^d* . 


For large numbers n,, using the asymptotic behavior K(k) ~~ — log Vl ~ k for k — *• 1, one obtains 

r. — \. v -k n -\ log n , (74) 

it 

which, a part from a logarithmic correction, has the same form of the preliminary result (53). 

5. NUMERICAL RESULTS 

The quantum evolution of input number eigenstates for the Hamiltonians (49-50) is evaluated 
numerically, taking advantage of the tridiagonal forms (28) and (40), which now read 

= \/n(j - 2n + 1)(j - 2n + 2) , (a = s„) , (75) 

- n + 1) , (a = s D ) ■ (76) 


C~L 


(75) 

(76) 





Figure 3: The maximum effective gain G w (corresponding to the conversion time T m of Fig. 2) 
for parametric Hamiltonians (49-50) (PNA on the left, PND on the right). The squares are for 
amplification /duplication, the circles for the inverse operating mode. 


The evolution of the output signal has been checked using the numerical results in Ref. [11]. In 
Fig. 1 we report a sample of the evolution for the PND. The time-dependence is periodic or 
nearly periodic for very low input photon numbers n, , whereas it becomes more and more irregular 
(essentially irreversible) for increasing n,. Qualitative differences between the direct and the inverse 
operating cases are evident. In the direct case the output signal exhibits maxima corresponding 
to high noise level, whereas low noise occurs only for depleted signal. In the inverse case, on 
the contrary, the first occurrence of a local maximum for the signal coincides with the absolute 
maximum, whereas the relative noise is always well below the subsequent values (this gap being 
an increasing function of n,). The conversion is never complete in both cases, however, it is more 
efficient in the inverse operating mode, due to the low noise at the output. The conversion time r. 
has been identified as the time corresponding to the first local maximum of the signal (in the direct 
operating mode this could be slightly lower than the absolute maximum). The same features in the 
time evolution can be found for the PNA approximated by the conversion Hamiltonian (49), with 
analogous differences between the direct and inverse operating modes. 

In Fig. 2 the conversion time r m is plotted against the input number n,, for both Hamiltonians 
(49) and (50). The direct and inverse operating modes lead to two different curves, the former 
corresponding to longer conversion times r. (a part from some features which are peculiar of the 
deamplifier for low inputs n,, and are reminiscent of the fractional behaviour (68)). The mean field 
approximation, which is pertinent only to the inverse operating mode, is reported for comparison. 
A good agreement is found for large n t , better for the PND than for the PNA. For large n, numerical 
best fits give power-low behaviours of the form r. ~ n" Q , with a .4 or smaller. 

In Fig. 3 the maximum effective gain G. (corresponding to the conversion time r. in Fig. 2) 
is reported. One can see that parametric conversion when used as a gain-two number amplifier 
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Figure 4: Effective gain G. for the deamplifier with G = 2: circles and full line describe the 
parametric conversion (Hamiltonian (49)); triangles and full line describe the intensity-saturating 
Hamiltonian (45); dot-dashed line corresponds to the ideal deamplifier (5). 

leads to an effective gain G m which is a decreasing function of the input signal n f , approaching the 
value G . ~ 1.28 for large n M well below the ideal gain. In a similar fashion the effective gain of the 
duplicator G . = (n a (r„))/(n c (0)} tends asymptotically to G* ^ .78. The inverse operating mode, 
on the contrary, behaves quite well, the deamplifier achieving the ideal G. = 1/2 gain and the 
recombiner G. = 1 in the large n, limit. The deampliiier gain is compared with the ideal one (5) in 
Fig. 4, where also the intensity-saturating case (45) is reported [notice that in the direct operating 
mode the intensity-saturating Hamiltonians (45-46) lead to ideal behaviour]. 

Finally, in Fig. 5 the output Fano factors F m at the conversion time r. are plotted. It is evident 
that parametric conversion lead to noisy PNA and PND, with F m ~ and exponent 0 slightly 
lower then 1: thisj:orresponds to an output signal-to-noise ratio whi ch is slowly (logarithmically) 
vanishing. The number deamplifier and recombiner are better approximated, with F m 2s .13 for 
large n t : on the other hand, the intensity-saturating Hamiltonian (45) leads to vanishing F. for 
large n, (F. is exactly zero for even n,.) 

6. CONCLUSIONS 

We end with some remarks on physical interpretation of numerical results. We have seen that 

parametric Hamiltonians (45-46) are not good candidates for number amplification/duplication 
devices, whereas they could be profitably used to achieve approximate number deamplifica- 

tion/recombination. Here we emphasize that the source of noise in the simulated number devices is 

the socalled spontaneous parametric emission. 6 * * * * * 12 As a matter of fact, as explained in Ref. [11], the 

Hamiltonians (45-46) are formally similar to the Hamiltonian of a laser amplifier: in particular, Eq. 

(46) can be put in correspondence with the Hamiltonian describing a cluster of N two-level atoms 
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Figure 5: The Fano factor F. at the conversion time r. of Fig. 2 for parametric Hamiltonians (49-50) 
(PNA on the left and PND on the right). The squares are for amplification/duplication, the circles 
for the inverse operating case. The triangles correspond to the intensity-saturating Hamiltonian 


interacting with one (resonant) mode of radiation 

Hi <x a'j- + aJ+ , (77) 

where J a = <r* are the collective spin-flip operators for atoms. In fact, the angular momentum 

operators can be represented in terms of the two mode-operators b and c as follows 

J + =bc\ J.= 6 f c , j, = i (c f c - *♦&) , 

(78) 

J = \(Jc + b'b), J-y. 

When operating as a PND the Hamiltonian (46) acts on input states with n„ = n*: in the direct 
operating mode one has (n a )o = 0 and (n c )o = n,, whereas in the inverse (» a )o = i n e) u = 0, 

namely \M\ = J in both case: this is exactly the spontaneous emission limit for the parametric 
converter (as opposed to the noisless coherent superradiant limit corresponding to M — 0). Thus, 
in conclusion, both the output noises and the nonideal effective gains are signs of the spontaneous 
parametric emission in the converter. 
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